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ABSTRACT 

We consider an ultra-relativistic wind consisting of electron-positron pairs and 

photons with the principal goal of finding the asymptotic Lorentz factor for zero 

baryon number. The wind is assumed to originate at radius rj where it has a Lorentz 

factor 7j and a temperature T, sufficiently high to maintain pair equilibrium. As r 

increases, T decreases and becomes less than the temperature corresponding to the 

electron mass me, after which non-equilibrium effects become important. Further out 

in the flow the optical depth r drops below one, but the pairs may still be accelerated 

by the photons until r falls below ~ 2 x 10~^7^'^^. Radiative transfer calculations 

show that only at this point do the radiation flux and pressure start to deviate 

significantly from their blackbody values. The acceleration of the pairs increases 7 by 

a factor ~ 45 as compared to its value at the photosphere; it is shown to approach 
7oo ~ 1.4 X 103(ri/10''cm)V4^3/4^,^^^_ 

The limit of zero baryon number is a good approximation when the 
mass injection rate M in the flow is below a critical value corresponding to 
(E/M)cfi ~ 5 X 10''(rj/10®cm)Tj/me for fixed energy injection rate E. For large 
baryon loading, E/M < {E/M)c,m ~ 350{ri/l0^cm)^/'^-ff^Ti/me, the asymptotic 
Lorentz factor is 700 ~ E/M. Surprisingly, increasing E/M from (E/M)c^m to 00 only 
increases 700 by a factor ~ {nip/me)^/^ ~ 6.5, less than an order of magnitude. 



Subject headings: gamma rays: bursts - hydrodynamics - radiative transfer 



1. Introduction 

The release of a large amount of radiative energy into a small volume can lead to the formation 
of a fireball, a dense fiuid of radiation and particles that expands under its own pressure. Fireball 
models have become the accepted framework for understanding gamma-ray bursts at cosmological 
distances and their afterglows (Paczyhski 1986, 1990; Shemi &; Piran 1990; Meszaros & Rees 
1993, Piran 1997). Paczyhski (1986) and Goodman (1986) originally considered the possibility 
that fireballs could originate in the collision of a pair of neutron stars in a binary star system 
coalescing as a result of gravitational radiation reaction. (See also Naryan, Paczyhski &: Piran 
1992.) In this picture, the thermal energy released in the collision, ~ 10^^ ergs, is radiated as a 
neutrino-anti-neutrino burst. A fraction of that energy may be transformed into electron-positron 
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pairs above the surface of the neutron star (Goodman, Dar &; Nussinov 1987). There are now 
additional proposals for the origin of fireballs (e.g. Paczyhski 1997; Fuller &; Shi 1997; Pen, Loeb 
& Turok 1997). 

Close to the radius at which energy is injected, the resulting wind is opaque. The radiation 
energy that is initially trapped can escape in two different ways further out in the flow: When 
the plasma becomes optically thin, radiation streams freely to the observer (Paczyhski 1990). 
Alternatively, if there is a significant baryon contamination in the fireball, it can become matter 
dominated before radiation escapes. The matter will increase the opacity and, more importantly, 
convert part of the radiation energy into bulk kinetic energy (Shemi k, Piran 1990). Interactions 
between the expanding atmosphere and the surrounding matter provide a way to convert the 
kinetic energy in the baryons back to radiation at the resulting shock front (Meszaros & Rees 
1993). Internal shocks due to variations in the velocity of matter is also proposed as a way to 
dissipate kinetic energy (Rees & Meszaros 1994). Fenimore (1997) recently pointed out that the 
observed temporal structure of gamma ray bursts severely constrains the proposed models of 
energy conversion by relativistic shocks; conceivably internal structure and shocks can account for 
some of the observed variability (e.g. Kobayashi, Piran & Sari 1997). 

In the fireball models invoked to explain the afterglows of GRB 970228 and ORB 970508, the 
bulk Lorentz factor of the outflow is 7 ~ 100 — 1000 before the outgoing shell is slowed significantly 
by sweeping up matter from ambient gas (Wijers, Rees & Meszaros 1997; Waxman 1997; Waxman, 
Kulkarni &; Frail 1998). From a theoretical point of view, values of 7 in this range yield acceptable 
estimates of burst duration and (with additional assumptions) characteristic photon energies. It is 
presumed that most of the energy originally in the fireball converts to kinetic energy long before 
deceleration begins, so the (baryon) rest mass of the flow must be nonzero: M = -E/7, where M 
and E are the rest mass and total energy injection rate of the flow. Precisely how such small but 
nonzero M /E arises is not clear yet; nor is it obvious whether 7 = E/M can be much larger or 
smaller than 100-1000, the values that seem necessary for modeling gamma ray bursts. 

In this paper we reconsider the original steady wind problem first solved by Paczyhski (1986) 
but for M = 0. At first sight, one might think that the result would be 7 — > 00. However, 
the failure of equilibrium at low temperatures (once the pair density falls suflficiently so that 
annihilation becomes slow) leads ultimately to a finite 7. As we shall see, in regions where the 
temperature is greater than the electron mass, the fireball is very optically thick because of 
the large number of electron— positron pairs. As the radius increases and the local temperature 
decreases, the deviation from equilibrium in the number density of pairs becomes significant. A 
little further out in the flow the optical depth falls below unity. However, the remaining pairs 
are still heated and accelerated considerably via their interactions with the radiation fleld; in the 
rest frame of the pairs, the radiation field itself remains close to the blackbody form long after 
the fireball becomes optically thin. The radiation spectrum detected by a stationary observer 
would not be blackbody, however, but also differs from the power-law spectra of GRBs. Although 
our calculations pertain to a steady, spherical wind, the general result that 7 is finite even at 
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zero baryon mass may be true as well for thin shells emitted from impulsive energy release. We 
consider non-spherical perturbations around our wind solutions in Section as we shall see, 
some memory of surface 'hot-spots' may persist out to the photosphere. 

The dynamics for winds with sufficiently small baryon number are similar to M = 
outflow. A larger baryon number will increase the inertia of the flow, and so we expect the 
radiative acceleration and therefore the asymptotic Lorentz factor 700 to be reduced. However, 
as the baryon number increases, the photospheric radius tends to increase as well, with a 
corresponding increase in the Lorentz factor at the photosphere. These two factors combined 
result in a surprisingly small variation in 700: For fixed E, as E/M decreases from 00 to 
{E/M)c,M ~ 350(ri/10^cm)^/^7^^^Tj/me, 700 is only reduced by a factor of ~ 10. For even smaller 
E/M, the baryons will dominate the energy of the flow even inside the photosphere, which results 
in a final Lorentz factor 700 ~ E/M. 

We present a detailed analytical model of the dynamics of the fireball in Section |^. This 
includes approximate results for the asymptotic value of the Lorentz factor and of the energy 
content in the pairs relative to that of the radiation, based on the initial temperature and initial 
velocity of the flow. Furthermore, we estimate where the pairs go out of equilibrium, the position 
of the photosphere, and the radius and optical depth at which the radiation fields start to deviate 
from their blackbody values. In Section |^ we show results from a numerical calculation to which 
the analytical model is compared. Next, based on the equation of radiative transfer, the comoving 
frame photon distribution function is shown to be very close to blackbody even out to quite 
small optical depths in these ultra-relativistic flows. The dynamical importance of baryons is 
described in Section |^ There we obtain approximate results in different regimes characterized by 
the amount of baryon loading, and we integrate the dynamical equations with baryons included 
in order to show how they affect the asymptotic Lorentz factor. Finally, we discuss qualitatively 
several possible extensions of our model, including the effects of (1) muon pairs and even nucleon 
pairs which could be present at sufficiently high temperatures; (2) a temperature anisotropy at 
the inner boundary, and (3) magnetic fields. 



2. Equations 

Consider a situation in which a large amount of energy is released into a compact region 
and the resulting relativistic outflow expands into vacuum. For simplicity we assume spherical 
symmetry and a stationary flow; gravity may also be neglected since we are interested in 
super-Eddington luminosities. (Perturbations away from exact spherical symmetry will be 
considered in Section ^^.) The cross-sections for absorption and scattering, cja and cr^, are taken 
to be constants; letting ctq and as cover a wide range of values in the numerical calculations 
facilitates a qualitative understanding of the effects of matter-radiation interactions. Throughout, 
units for which ks = h = c = 1 are used. 
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The temperature at the radius where energy is released is assumed to be high enough for 
pair creation to occur, and the flow is very optically thick close to its inner boundary. When 
the optical depth r is 3> 1, the radiation field in the comoving frame will be close to blackbody. 
In a static or slowly expanding atmosphere one would expect the deviations from blackbody in 
the radiation field to be of order 1/r. However, in an ultra-relativistic (Lorentz factor 7^1) 
expanding atmosphere, the corrections to blackbody actually vanish to first order in 1/r and 
lowest order in I/7. This will be discussed further in Section 0. 

In the opposite limit, the optical depth approaches zero and the radiation streams freely. In 
this case a photon preserves its frequency in the lab frame as well as the quantity (1 — ^^)r^, 
where /i is the direction cosine relative to the local radial direction. As the radial coordinate 
r gets large, the intensity therefore becomes increasingly sharply peaked in the outward radial 
direction; this is a purely geometrical effect (Hummer & Rybicki 1971). In a static or slowly 
moving atmosphere one would expect the radiation quantities to develop large deviations from 
their optically thick values soon after reaching the regime where r < 1. However, we shall find that 
in a rapidly expanding atmosphere the radiation field maintains its equilibrium form long after 
it decouples from the matter. This is similar to what happens in expanding universe cosmology, 
where the background radiation preserves its blackbody spectrum after decoupling. Specifically, 



as will be discussed in detail later (Section 3^), we find that the radiation fields start deviating 



from blackbody radiation in the comoving frame only when r < 2 x 10 ^7^''^. 

This suggests a simple model for the radiation quantities: For r > 1, the radiation field is 
approximated as a blackbody distribution in the rest frame of the flow; for r < 1, the radiation 
energy, flux and pressure are given by the free-streaming approximation. We also need equations 
for the matter, characterized by its temperature T, Lorentz factor 7, and pair number density 
Ue- Energy and momentum conservation give two dynamical equations, whereas the Boltzmann 
equation determines the number density for the electrons and positrons. 



2.1. Equations for r > 1 

We assume that the flow originates at radius rj where it has a Lorentz factor 7^ and 
temperature Tj. The basic equations describing the flow are found from the energy-momentum 
tensor, 



I pV/^. (1) 



(27r)3 

p" is the four-momentum, / the distribution function, and E = \J-p^ -|- with m being the mass 
of the particle contributing to T"^. For the matter (e^ pairs for most of this paper) 

T'S = (PM + PmWU" + PMg'^f', (2) 

where pM and Pm are the energy density and pressure measured in the rest frame of the flow 
and is the matter four- velocity. We also need the following expressions for number density n, 
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energy density p and pressure P: 



" = / f{p,T)Wdp; (3) 



2 

(2V) 
2 

(2V)S 



= I \lv^ + nn?f{p,T)4.T^p'^dp and (4) 



P 



The non-zero components of the energy-momentum tensor for the radiation are 
TOO = ^0,7^ = FcT^i = Po and T^^ = T^^ = (l/2)(Po - Po) in the comoving frame, 
where the 1-direction is along the flow velocity; Lorentz transforming to an arbitrary reference 
frame gives 

T^o = ^''[Eo + 2vFo + v^Pg] (6) 
= ^^[v{E, + P,) + {l + v^)Fo\ = T'^ (7) 
T}1 = -i\Po + 2vFo + v^Eo] (8) 

Tf = l{Eo-Po)=T^^. (9) 

Below, we call Eo,Po and Fq the (comoving frame) radiation energy density, pressure and flux 
(e.g. Mihalas & Mihalas 1984). 

The equations describing the flow are given by T.^^ = 0, where T"'^ = + T"^ . With the 
assumptions of spherical symmetry and steady state, one gets 



j,2 fj^j. 



{PM + Pm + Eo + Po)7V7^ - 1 + (27 - l)i^c 



(10) 



and 



-MA{pm + Pm + Eo + Po)(7^ - 1) + 27y7^Po]}+^ + ^ + = 0. (11) 

ar * dr dr r 

For large r, the radiation has a blackbody distribution in the rest frame of the flow, and 
Eo = TT^r^/lS, Po = and Pq = Po/3. 

When the temperature in the rest frame of the fluid falls below nie, the equilibrium 
+ e~ ^ 27 can no longer be maintained. Consequently, the number densities of electrons 
and positrons will deviate from their equilibrium values and must be found from the Boltzmann 
equation, 

= (^) ■ (12) 

TTlg dx"^ V dt j collisions 

Integrating over ^p we get, in spherical symmetry. 



collisions 
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For the process 



7 + 7 , 



(14) 



the cohision term is obtained by integrating the relevant matrix element times the factor 
{f-yfy — f-f+) over phase space. Deviations from pair equilibrium become significant when T/rrie 
is small, so we may approximate /j <C 1 (i.e. fi = exp [— {Ei — nt) /T]). We can therefore neglect 
Fermi suppression of electron and positron final states and stimulated emission of photons in the 
regions where equilibrium fails and the pair density 'freezes out'. As will be seen later (Section 
3.2), equilibrium is maintained until T is ~ O.OSmg. With these approximations. 



due 
dt 



collisions 



n 



e,eq 



(15) 



where eq is the equilibrium number density of electrons; using this result, equation (O) can be 



written 



1 d 

J.2 fj^j. 



n 



e,eq 



Svensson (1982) constructed a useful approximate expression for (cTann^) 

„2 



vrr; 



1 + 2(r/me)2 / In [2r]ET/me + 1.3] 



(16) 

valid at all temperatures: 
(17) 



Here is the classical electron radius, rjE = exp(— C^;) and Ce ~ 0.5772 is Euler's constant. For 
the purpose of our analytical estimates we only need to know that {(Ta,nav) ~ vrrg for T <C 



2.2. Equations for r < 1 

For small r, a good approximation to the radiation fields is found by assuming that photons 
stream freely from the photosphere to the observer. However, as the optical depth drops below 
one, energy and momentum may still be deposited by photons in the pair wind. Even though the 
energy and momentum deposition may be small compared to the energy-momentum content of 
the escaping photons, it may be large compared to that of the pairs. The dynamical equations 
must therefore take this interaction into account. We shall see that substantial acceleration of the 
pairs occurs at r < 1. 

The equation of radiative transfer can be written schematically as 

k^^f^ik,x) = C^. (18) 

f^{k,x) is the phase space distribution function of the photons, and the Lorentz invariant collision 
terms C-y result from interactions with the electrons and positrons in the fluid. 

We are interested in three processes: — 7 scattering (e^ + 7 — > + 7) , bremsstrahlung 
and its inverse (e^ + ^ + + 7), and pair annihilation and creation (e"^ + e~ ^ 27). We 
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use a phenomenological model for the collision terms, evaluated in the rest frame of the flow. For 
scattering the collision term is 



f-y{k,x) + d h' g{n' h)f^{k',x) 



(19) 



where is the photon energy in the rest frame of the flow, and we assumed that the scattering is 
elastic (A;"^ = A;"). The probability distribution for direction changes in the rest frame is g{h' h). 
The cross section for scattering, as, is taken to be a constant. The collision term for absorption 
and emission is 

Ca = 2neaa~k'>[-f^{k,x) + fBBikVT)], (20) 

where fsBik^ /T) is the equilibrium photon distribution function at temperature T. The absorption 
cross section da is also assumed to be constant. 



k'^C.y. (21) 



Now take the first moment of the Boltzmann equation (|18|): 

d^k 



{2TTf J kO 



d 



(27r)3y fcO 



We recognize the left hand side as T"^!^- The right hand side is the radiation force, which we 
denote G". The physical interpretation of the equations T^!^ = (or T'^.^ = — G") is that —G^ 
equals the net rate of radiative energy density input into the matter, while — G* is the net rate of 
radiative momentum input (Mihalas & Mihalas 1984). For the scattering term, taking the first 
moment over Cs as indicated in equation (|2^) gives 

^8,0 = G's.O = G's,0 = 0; = -2ne(TsFo, (22) 

where lower index denotes that these quantities are evaluated in the comoving frame. The 
components of the radiation four-force corresponding to absorption and emission are 

G° = -2n,aa[E^ - U,^{T)] G^^ = -2n,aaFo G^ q = G^ ^ = 0, (23) 

where Ueq{T) = Finally, for pair annihilation and creation we use the interpretations of 

G" to write down directly 

G% = 2me(aa„n^^7)[n^ -<eq(r)] G^o = [i = 1,2,3]. (24) 

The energy loss of the plasma per annihilation is (7_ + 7+)me. Svensson (1982) gives an 
approximate expression for the pair annihilation cooling rate: 

^'""""''^^ " 1/(1 + 6T/me) + (T/me)/ {In [2r]ET/m, + 1] + 1/4} ' ^^^^ 
In the non-relativistic limit the energy loss is dominated by the rest mass energy, and one gets 

(0-annW7) ~ TTr^ for T < TTig. 



-8- 



Using the energy momentum tensor for the pairs and the components of the four force 
G"' as given above, the dynamical equations T^,^ = — G° reduce to 

[PM + Pm)-t] \- -T, — = Gq (26) 

a in r a m r 7 

and 

^ , , p„, , , p„) (i - ^) ^ - -«87^- (^^) 

Up to this point, no restriction to small optical depths has been invoked, and the dynamical 
equations can be applied at any r, provided f^{k,x) is determined. For r < 1, wc assume 
that radiation streams outward freely. This is only marginally accurate at r ~ 1, but becomes 
progressively more precise with decreasing r. However, since the photon intensity is sharply 
forward-peaked in the lab frame, this approximation is better than one might have expected even 
at T < 1. A non-interacting photon preserves its frequency in the lab frame as well as the quantity 
(1 — fP)r^, where fi is the direction cosine relative to the local radial direction. The phase space 
density fj{k,x) is conserved along any photon ray path. Consequently, we need only to relate the 
photon frequency and direction cosine in the comoving frame at different points in the flow to find 
Eq, Pq, and Fq at small r. In the extreme relativistic limit, photons are highly beamed in the 
forward direction as seen in the lab frame. It is therefore convenient to define the quantity 

77 = 27^(1 -/x), (28) 

in terms of which the direction cosine in the rest frame is approximately /ij-f ~ (1 — r7)/(l -|- 77). As 
T] spans the range [0, cxd], ^Urf spans the range [—1,-1-1]. Prom 2(1 — yu)r^ = •qr'^/j'^ ^ constant, we 

i+M:;'-(i-M:,")c(r) 
'^-i+m;;'+(i-m;;')cm ' ' 

where 

C(r)=(^)'. (30) 

The label '1' denotes a reference point in the flow, taken to be r ~ 1 below, where we shall assume 
that the distribution function is known and beyond which photons stream freely. Note that if 
C(r) = 1 then /Xj-f = /^rf ' which means that the direction cosine relative to the local radial direction 
does not change as long as 7 oc r; however, as ({r) 0, the distribution becomes increasingly 
peaked around /Xrf = 1- We shall see that 7 oc r is maintained to quite small values of r. 

In a similar way, the photon frequency k^i in the rest frame of the flow is 

fc(l + ^) ^ , (1) lif l + V 
27 7 U + m 

or, introducing /irf rather than 77, 



k.^ = -^^ = Kt^-{T^]^ (31) 



kri — k^t ^ 

7 



l + M;j^ + (l-M?)C(r 



(32) 
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Now assume that the photon distribution function is f-y{k^]\ at ri. If the photons flow 
freely outside ri, then at any larger value of r, the Liouville theorem implies that the phase space 
distribution is still f^yik^^K ^J-^^^)- Only the mapping from ri to r is needed to find Eq,Pq and Fq. 

Suppose the distribution function for the radiation at ri can be expanded in terms of 
orthonormal Legendre polynomials: 

= MkS') + ,^Mki^') + fM^) (33) 

Defining 

h= j dqq^foiq) 2Ji = J dqq^h{q) and lOKi = J dqq^f2{q), (34) 
the radiation energy, flux, and pressure can be expressed as: 

^0 = y (^)'{/i[l+C(r)+C'(r)]+Ji[l-C'(r)]+i^i[l-CF} (35) 

Fo = f i^^y {h[l - CHr)] + Ml + C\r)] + K,[l - e]} (36) 

Po = fi^^y{h[l-ar)+e{r)] + Ml-CHr)]+K,[l + C?}. (37) 

We can check that in the stationary frame r'^F{r) = constant, consistent with energy conservation. 

For our purposes, the reference radius ri will be taken to be the position of the photosphere. 
The radiation field is then given by an isotropic blackbody distribution for r < ri . (In Section we 
estimate deviations from a blackbody.) We therefore specialize to /i = /2 = (i.e. Ji = Ki = 0) 
in the following discussion. 

As long as = Fq = 3Pq and Fq = 0. However, as Ci^) ^ 0, Pq ^ Eq and Fq — > Eq; only 
when 7 — > constant do Pq, Eq and Fq decrease like r~^. Note that one can rewrite the prefactor in 
equations (35-|37|) as 

vw) -^(^' ^''^ 

consequently, when Qir) = 1, the rest frame energy density and pressure both decrease like r~^. 
This is so even though there is no interaction with matter. The situation is the same as in 
expanding universe cosmology, where the radiation field maintains its equilibrium form after it 
decouples from the matter entirely. We see from the above that when the optical depth is small, 
the critical function that determines the approach to Fq/Eq = 1 and Pq/Eq = 1 is C(^)- 

For the pairs, we always solve the annihilation equation derived in the previous subsection. 
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3. Analytical Model 

In this section we will do back-of-the-envelope calculations and obtain approximate results 
and scaling laws in four different regimes, each characterized by their local temperature T and 
optical depth r: First, in the high temperature limit (T S> rUe), electrons, positrons and photons 
are in equilibrium, the optical depth is very large and we get simple power laws for 7(r) and 
T(r). As the temperature drops below mg, the number density of electrons and positrons starts 
to deviate from its equilibrium value and will be calculated approximatively from the annihilation 
equation. In the region where the optical depth r < 1, we estimate how the radiation fields affect 
the velocity and the temperature of the matter. Finally, we calculate the asymptotic values for 
7, Tr and ngr^ for r ^ 0. A summary of the scaling laws can be found in Fig. |l[ The validity 
of the simple analytical model presented here will be checked against a numerical solution of the 
equations in the following section. There we will explore models with a range of different initial 
temperatures and Lorentz factors. 



3.1. T > 1 and T » m-e 

When T ^ mg, the radiation energy Eq is given by the equilibrium Planck energy density, 
Eq = (7rVl5)T^ = aT^ = SPq. The radiation flux is negligible since the optical depth r ^> 1; we 
will check the consistency of this assumption later. (See Section 3^ for calculations of the optical 
depth.) 

Introduce the notation 

p = PM + Eo = Pe~ + Pe+ + Eo and P = Pm + Pq ^ Pe- + Pe+ + Po- (39) 
Energy conservation, equation (10), now simplifies to 

d_ 

dr 



r\p + P)-id-f'^ -l\ i.e. L = 47rr^7J72 - 1(/) + P) = constant. (40) 



The Euler equation, equation (|llD, becomes 



(iln7 dP ,, , 



where we made use of energy conservation (equation ||40|| ). 

When p and P depend on T alone, dP/dT = (p + P)/T (e.g. Weinberg 1972, Section 15.6), 
and the fiow equations give 

^^"/^^^ =0; -fT = constant = -/iTi, (42) 
amr 

and 

dln7 ^ 2 ^ 2(7^ - 1) 

dlnr (rdp/(ir)/(p + P)- 72/(72-1) 272-3' ^ ^ 
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The last step in the above equation is vahd for p = 3P oc T^. For 7 1, dlwy /dlnr = 1, i.e. 

7 tx r; T oc -; 7T = 7jTj for T » mg (44) 
r 

(e.g. Goodman 1986; Paczyiiski 1986). Note that equation (^) requires a minimum initial velocity 
for the flow to lift off: Demanding dlwy /dlnr > implies that 7^ > ^3/2. 



3.2. T > 1 and T nie 

When T ^ rrie, Pm{= S-Pm) is comparable to -E'o(= 3Po) and also proportional to T^, but for 
T ^ nie, pm ^ Eq. In this subsection we want to explore what happens when T < nie- Define 

In general, pM and P/vf should be found from equations (Q) and (^). 

We can now calculate the radius where the temperature equals the electron mass, using 
conservation of energy: 



ri 



. h{T,) _ ( iV'^m, / h{T,) (._]_\'^T, 

\ ImVl^Kme) \ 1! V him,) ^ I 72 ] m. 



Here 7^ ^ 1 and 7m"^e ~ 7j7i was used. Note that the pair density is still close to equilibrium 
when T ~ rup, which will be verified later in this subsection. 



as 



As T drops below rrie, the equilibrium density of electrons and positrons falls off exponentially 

:(meT)3/2exp(-me/r). (47) 



Consequently, for T <C me, Pe,Pe Eq,Pq, and p ^ Eq 3P ^ SPq « (7r2/15)T^. But then 
we get the same scaling laws as in the high temperature limit (see the justification of equations 
©-HI): 

1 

7 (X r; T oc — ; 7T = constant for T <C mg. (48) 
r 

As long as equilibrium is maintained, the equation of state of the flow is purely a function of 
temperature, and the Euler equation implies that jT = constant. At the same time, energy 
conservation requires 



L PS 47rr2/i(T) 7^/72 - 1 ps Airr'^i-fiTifT'^a = constant. (49) 
Thus, r^T^a = constant, resulting in the scaling laws 

7 ss 'Jiir/ri) and T w Ti{ri/r) for T > me] (50) 



7 ss 7i(r/ri)y47ll and T Ti{ri/r)y/lT/l for T < mg. (51) 
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The change in slope is due to the fact that pairs contribute httle to p + P once they become 
non-relativistic. This remains true for smah T < mg even though are far more plentiful than in 
pair equilibrium. 

At somewhat larger radius than rm, deviations from pair annihilation equilibrium will become 
important. (Lee & Weinberg (1977) considered the analogous problem for heavy leptons in the 
early universe.) In order to estimate where this happens, let us assume that the deviations are 
small, i.e. 6ne = rig 
to get 



ne,eq f^e.eq- Then expand n; 



n 



e,eq 



Irie^eq^rie in the pair equation (|l^) 



n, 



e,eq 



r W7. 



Since 7 oc r as long as equilibrium is maintained, we get 



7 



d In Up p„ 

2 _j_ *^5^M 



din r 



(52) 



(53) 



From the exponential form of rie^eq, dlnne^eq/c^lnr = —dlnne^eq/dlnT « — m-e/T for T <^ nie. 

Equilibrium fails when 5ne becomes comparable to rie^eq- Let r^q (and T^q) be the radius (and 
temperature) where Srie = Ue^eq] then 



(27r)3/2 



[m 



eqj 



Icq rUe 
'^'^eq'^^ e '^eq 



(54) 



to leading order in mg/Teq. Relate Teq and r^q by TeqVeq ~ y/TT/4merm, and use 
7eq/?^eq ~ \/4/ll7j/ri (from the limiting cases of h{T), assuming Teq <C rUe and Tj > rUe 
and 7-^/72 — 1 PS 7^ in the energy conservation equation). We then find 



(27r)3/2 



11 \ 3/4 

4 



3/2 



' eq 



which has the approximate solution 



Teq/rm ~ -y/ll/41n 



exp(- j4/llreq/r^ 



/2/7r(ll/4)V4(,2,^^3)/^^ 



(4/ll)(7i/r,)r, 



eq 



27rr2 



11/4(5/2) In 



' eq 



(55) 



(56) 



Note that req/r^ does not depend on the initial temperature Tj and is only weakly dependent on 
7j. For Tj = 10^ cm, the numerical value is req/r^ ~ 33 — ^ll/41n7j. The temperature at r^q is 
approximately Teq ^ yJll/Amerm/feq ~ 0.052me. 



3.3. T ~ 1 and T <^ nie 

At Teq the energy density of the flow is dominated by the photons: Since Teq n^e have 

1.1 X 10"Si- (57) 



^0,eq (vrVl5)r4 



eq 
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As long as the flow remains photon dominated and the radiation field stays close to its blackbody 
form, 7 oc r. 

Now, let us estimate the optical depth in the fiow. In the rest frame of the flow, the density 
of electrons and positrons is 2ne = n_ + n+ and the cross section for interaction is {as + Ca)- A 
reasonable choice for the optical depth of the flow is the radial depth (Abramowicz, Novikov & 
Paczyhski 1991) 

— = 2ne{(Ts + (Ta)\ — ^ , (58) 

dr V 1 + 7 

where the last approximation holds for an extremely relativistic flow. We then find 

3 Vo-T or/ Jr l{r) 

As long as the radiation field remains close to a blackbody, ^{r) ~ ^eq(j'/i^eq)- Since > ne,eq 
outside Teq, we can approximate the pair equation as 

In other words, although pair annihilation continues outside req, pair creation becomes unimportant 
as the number of photons energetic enough to produce pairs decreases exponentially. The 
solution to this equation is 

ne{r)^ ^ "^^'""^ r ^f^V (61) 

l + (l/3)^47TT(req/r^) l-(req/r)^'' ^ 



for r > req and 7 oc r. Since ne(req) ~ ne,eq(req) + SUeireq) ~ 2(5ne(req) 
(4/ll)(7,/ri)(req/r^)/7rr2 



.2 



/X Stt 2 / 0-a . 0-5 \ req . . 3 /■°° dr 
r(r) — r^ 1 — -rielreajr^^ / — p r 

3 Kar cttJi^ ""Jr r^ {l + (l/3)v/4/ll(req/r^) [l - (req/r)^] } 



\3 



8 / ^\ ^ r _ (l/3)v/47TT(req/r^)(req/r)^ . 
3Ut cjtJ \ l + (l/3)v/47TT(req/r„) J ^ ^ 



at r > req. The photospheric radius rph may be estimated from r(r = rph) = 1, implying 

-1/3 

1 — exp 



' cq I \ ' ^ ^^'i 



3aT 



8 [(7a + as 



(63) 



For example, if Ua + o"^ ~ ax then rph/rcq ~ 1.4. This implies that the temperature 
at the photosphere is Tph/irie ~ (req/rph)req/me ~ 0.037. The optical depth at req is 

Teq ~ 5.4((7„ + as)/aT. 

Outside rph the photons stream nearly freely out to the observer. The escaping photons form 
a high energy gamma ray continuum since T^^ the temperature seen by a distant observer in the 
lab frame, is approximately constant in the flow out to rph. 



- 14 - 



3.4. T < 1 and T <. irie 

For r < 1, energy and momentum may still be deposited by photons in the pair wind. Thus, 
the pairs can be accelerated and heated via their interactions with escaping photons even though 
the wind is not opaque. Because the energy and momentum in pairs is smaller than in radiation 
at r ~ 1, even small depositions of energy and momentum may affect the pairs significantly. In 
this regime, each electron or positron experiences a force 

/;f ~ {cTa + <rs)Fo (64) 

in the flow rest frame. (Since T <C rrie, there is little difference between rest frames of flow and 
individual electrons.) Transforming to the observer's frame implies 

d'y 

-r ~ (o-a + (Ts)Fo/me, (65) 
ar 

where we used dr = dt/^ and dt = vdr ~ dr. This can also be seen from the dynamical equation 
(p^. The appropriate flux in this regime is given by the free-streaming limit, 

fo = ^?f^)^.[l-CV)l. (66) 



3 \ vy 

assuming, for simplicity, that the radiation field at ri is isotropic, as will be the case when 
matching onto blackbody radiation at r ~ 1. Define z = r/ri and 

/ \ 4 

47r /fJa , CTs \ , 3 /i Zi f <^a . (ys\ U 



A^^I-+l±UaTrmt)— ,-^3.1x10' [- + - --^1^ . (67) 
Since 7/71 = (2;/2;i)v^, the equation of motion becomes 

At the photosphere, z = zi and C = 1- Since A ^> 1, will remain close to one until 2A{zi/z)'^ ~ 1. 
Thus, 7 oc r until z = z-y, where 

^ ~ AV4; (69) 

^1 

note that for zi = rph/rj, r^/rj ~ 2 x lO'^Tj/me for and (da + as)/c7T of order unity. In this 
simplified model, 7 700, the asymptotic Lorentz factor, for z > z^. Moreover ^00/ z.-^ ~ 7i/-2i> so 

7oo ~ 7iA'/^ (70) 

The resulting 'boost' 700/71 may be large: For example if 7j ~ 2 and aa + (Tg ~ cti then 
A ~ 4 X 10^, and the asymptotic Lorentz factor of the flow increases by a factor of ~ 45 as 
compared to its value at the photosphere! Specifically, by using previous estimates for the flow 
variables at r > 1, we get 

7oo ~ 1.4 X 10^7/ — + — T7^6 • Ci^^ 

me \(Tt (Tt J V10°cm/ 
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We are now in a position to justify our assertion that the radiation field wih stay close to 
blackbody in the comoving frame long after the optical depth drops below unity. In the free 
streaming limit the radiation energy, flux and pressure remain close to their equilibrium values as 
long as C — 1 is small, which will be the case until r ~ r^. For r > r^,7 ^ 700, and the right-hand 
side in the pair equation can be neglected. This is because annihilation becomes unimportant for 
r > because of the low number density of pairs, as can be easily checked. Thus, for r > r^, 
Uer'^ ~ constant and the optical depth at becomes T(r^) = w ne{r.y){r^ / 'yao){o'a + Cs). 
Using our previous results, we find that little deviation from blackbody radiation develops until 
r, ~ 1.7 X 10-5 (f^ + f^)'/^f^ 

Finally, we need to estimate the asymptotic number density of pairs. Equation (^Tj) gives the 
approximate number density for r^q < r < r^, where 7 oc r: 



eq ^ ' ^'71 

eq 



(4/ll)(7,/ri)(req/r™)(l/7rr2) (r, 



for r > Teq; 7 oc r. (72) 



1 + (l/3)y4/ll(req/r™) [l - (req/r)^ 
For z> the pair equation reduces to 

ne{r)z^ « n^(r^)zl^ ne.{r^)z^-^b^l'^ . (73) 
Using the approximate results for ne(r^), rm/fi, r^qlrm and rph/req, we find 

ne ~ 8.6 X 10^9 { ^) ' (7.)^/^ f^^^) ( L\ " em-3 (74) 
\me) \OT ot) \ri) 

for r > r^. 

The acceleration of the pair wind increases the energy outflow in pairs to a fraction 
Le/L PS 45r^7^pAf /ll7r^r?7?Tj'^ of the total luminosity. Using pu ~ 2mene and earlier estimates 
for the other quantities in the above ratio, we find 

LJL ~ 8.5 X 10-6 (7,)3/4 + ^\ (75) 

VcTT fr/ 

for r > Try We conclude that for reasonable 7^ the flow never becomes matter dominated. Almost 
all of the original L is carried away by photons. Notice that although the asymptotic pair density 
and Lorentz factor depend on both 7j and Tj, the value of L(,/L at r ^ depends only on 7j. 

To find the asymptotic behavior of the temperature of pairs, consider the entropy per particle, 
Tds = PMd{l/2ne) + d{pM /^rie). For T <^ rrie, Pm ~ 2neT and pM 2mene + 3neT. Since 
rig oc for large r we have 

Tds ^ -Tdr + -dr. (76) 
r 2 
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Adiabatic cooling, s=constant, would result in T oc r~^l'^ (as could also be obtained from T oc n^/^ 
and n oc r~^). However, energy is still deposited by photons in the flow even at r > r^; the 
amount of radiative energy input per particle in a time t,- ~ r /7 is 

dq = -Gltr/2ne ^ -GE]r/2ne7. (77) 

For r > r^, £^0 ~ -^0 ~ -Po and energy conservation implies that L ~ Attt'^j'^IEq + Pq + 2Fq] « 
167rr^7^£'o. Consequently, —Gq « 2ne(TaEo « 2neCraL/167rr^7^. Since Tds = dq, equation ( [76| ) 
takes the form 

dlnT 4 (TflL 

dlnr ^ ~3 2ATTj^rT ^ ' 

with solution 

T ^ + ^ (79) 

There are two cases to consider, depending on the relative sizes of the two 
terms on the right hand side in equation (|7^); their ratio is aaL/327r'j^rT ~ 
7.5 X 10^((Ta/(TT)7~"'^^^[((7a + ffg) / aT]~^^'^{ri / 10^ cui) . If this ratio is larger than unity, 
T oc r~^. If the ratio is much smaller than unity, T oc r"^/^ until a radius 
r/r^ ~ (87r7^r^r^/aaL)3 ~ 3.3 x 10-^HaJaT)-^jf^[{aa + as)/aT]^/^{ri/10^cmy\ after 
which T oc r~^. This result could also have been obtained from equation (p7|). 



3.5. Positronium formation? 

To this point, we have considered annihilation in flight only. At large radii, pairs 
cool considerably, and it is possible that recombination to positronium occurs, leading ultimately 
to destruction of the pairs. Like hydrogen recombination, positronium formation would proceed 
via a free-bound transition to an excited state, followed by a radiative cascade to lower energy 
states, with increasing probability of pair annihilation. However, in order for this process to be 
important, the rate at which positronium is formed must exceed the expansion rate of the flow. 
To estimate the rate of positronium formation, we use the results from hydrogen recombination 
with the reduced mass ~ mg replaced by me/2. 

It is straightforward to show that the 'gross' recombination rate - ignoring radiative 
ionization, which leads to a lower 'net' recombination rate - is small compared to the expansion 
rate Fexp ~ j/r at all points outside r^q. Recall that at req, the electron temperature is 
Teq ^ 0.052me ^ = Xpos, the ionization potential of positronium. As long as T ^ Xpos, 

the recombination coefficient for the nth excited state of positronium is r„ = ne{crv)n, where (e.g. 
Rybicki & Lightman 1979) 

, , 32^^ a^r^ /me f 2Tn^\ 
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At Teq, we find Ti/F exp ^ 10 "^'j for Tj > 1, the ratio Tn/^exp is even smaller. Moreover, 
Tn/rexp oc UeT^^/'^r /j OC r~'^/^ outside req, but still at radii where 7ocr, Tocl/r and T > Xpos- 
Thus, recombination is very slow in this regime. 

Ultimately, T drops below Xpos and the recombination rate (summed over all states) is 
approximately 

Tpos « 2.08 X 10"^3neA^/2(0.429 + (1/2) In A + O.mX-^^^) (81) 

where A = 1.33 x IQ-^nie/T and Tie is in cm ^ (e.g. Rybicki & Lightman 1979). In this regime, 
Fpos/rexp oc Uf.r/jT^^'^ oc l/r7^r^/^, since n^r'^j constant; thus, Tpos/r^xp continues to decrease 
with radius far out in the flow, at first oc r~^/^, and ultimately oc r~^/^. Thus we conclude that the 
formation of positronium is always very slow, and does not affect the pair density substantially. 



4. Results 

A large number of models were calculated, covering a wide range of initial temperatures 
(1 ^ Ti/nie < 1000), initial Lorentz factors (1.25 < 7i < 10), and absorption cross-sections 
(10^^ < OajoT < 1). The scattering cross-section was held fixed at (Ts/(Jt = 1) and the initial 
radius ri = lO^cm. Due to the stiffness in the pair equation, and in the dynamical equations for 
r < 1, we used a semi-implicit scheme for solving the system of equations. 

Figs § and |3| show the ratios rm/Ti, req/rm, '^ph/''^eq and r^/rph as functions of the initial 
temperature Tj and Lorentz factor 7j. The agreement between the numerical results (solid lines) 
and the approximations from the analytical model (dotted lines) is very good. 

Fig. ^ confirms the high asymptotic Lorentz factors as well as their dependence on initial 
temperature and Lorentz factor, 700 oc {Ti/me)''y^^^ ■ 



Equation ( |79[ ) implies that T oc 1/r for large r; this is seen in Fig. ^. The behavior of 

;/3 



(ilnT/dlnr can be understood from the dynamical equations T'^f = —G". For T <C rrie, pm ^ Pm 



and the energy equation can be approximated as 

3dlnT dln-y z ^ rUe Oa 

oTTI + ^ + ~J^ ^ -(f^Tr-^mJ — — 

2 otmr dmr 7 1 gt 



(82) 



The left hand side of this equation contains terms of order unity. If (Ta/'^T is small enough to 
make the right hand side much less than 1, we must have dhiT / dlnr ~ —4/3 — (2/3)dln7/dlnr. 
For r < r^, where 7 oc r, dlnT/dlnr will therefore tend to approach —2, as can be seen from Fig. 
I when aa/(JT = 10"*^. When r > r-y, 7 constant and dlnT/dlnr approaches —4/3. This is 
also seen in the graph, before heating from the radiation eventually forces T to be proportional 
to 1/r. On the other hand, if (Ja/cTT is large enough to make {z / ^){aTrim^){me/T){aa/ (Tt) 3> 1, 
Eq/ui'^ — Ueq/m^ must be very close to 0. Prom the free-streaming form of £'0 this means that 

oc r~^7~^. Consequently, for r < we have dlnT/dlnr ~ —1, while for r > r^, 7 ~ 700 and 
T oc r~^/^. This is seen from the plot of dlnT/dlnr for aa = (Jt- 
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Next, Fig. I serves as a check of our estimates for the asymptotic number density of pairs, 
and for the ratio of energy in pairs to the total energy. Again, the analytical model (dotted lines) 
agrees very well with the numerical results (solid lines). 

When r < r^, i.e. r > r^, the Lorentz factor approaches a constant value, as seen in Fig. |^. 
The dip in din ^/ din r for large optical depths happens when T ~ rrie, after which pairs contribute 
little to the total energy. This point is the transition between 7 ~ 7i(r/rj), when T ^ rrie, and the 
low-temperature limit 7 'yi{r / ri) 4/ 11 (see equations and |l5l| ). 



5. Robustness of Blackbody 

Here, we show that the deviations from blackbody in the photon distribution function are 
small. In the 7^1 limit, the corresponding spectrum observed in the lab frame has a slightly 
broader peak and a shallower slope for low frequencies compared to a blackbody. We also consider 
non-radial perturbations to the flow and resulting effects on the observed radiation spectrum. 



5.1. Deviations from blackbody in the comoving frame 



In this subsection we will give a more careful justification for the simple model we used for 
the radiation field. When the optical depth is large, the radiation is trapped effectively. Because 
thermal emission and scattering are isotropic in the rest frame of the flow, one expects the 
radiation field to be close to isotropic when the optical depth is larger than one. The deviation 
from isotropy can be treated as a small perturbation on the radiation field in this limit. From 



Section 2.2, the equation of radiative transfer is 



d 



7- 



The Lorentz invariant collision terms are 



fJk,x)+ / d'^n g{h' ^ h)fJk',x) 



and 



(83) 



(84) 



Ca = 2neaak'^[-f^{k,x) + fEBik^T)] (85) 

for scattering and absorption/emission, respectively. Below, we neglect pair production and 
annihilation since these processes have a negligible effect on the escaping radiation. 



Perturbing away from blackbody radiation, we write 

Uk,x)=hBCkVT) + 6Uk,x). 



(86) 



It is convenient to introduce the four- vector = k^P^'^ /k^ , where the projection tensor is 
P^^ = g^^ + U^U'^. This four- vector reduces to the unit vector along the direction of propagation 
in the rest frame of the flow. In terms of and k^ = —k^lf^, we have /c^ = k^in^^ + U^). 
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Now expand the correction to blackbody, 5/^, in Legendre polynomials: 



5f^{k,x)=A + B^nf' + C^, 



V 2 



(87) 



Introduce an explicit form for the scattering kernel, 

1 + aP2{n ■ n') 



g[n n) 



where P2{h,h') = [3{h ■ n')^ — l]/2. We will specialize eventually to Thomson scattering for which 
a = 1/2. 



With the expansions in equations ( pq ) and (p7[), the right-hand side of the transfer equation 
becomes 



n 



+ (Ts 1 - - 



+ Cp. 



(89) 



On the left hand side one gets 



(90) 

with ip{z) = zf^^{z). Equating equations (p9|) and (pO|), the results for A, and C^y can be read 
off. This allows us to find the energy momentum tensor for radiation, 



R 



The results are 
5^0 



d^^k 



k^k" 



M-k^Ux/T{x)) + A + B^n^ + C^a 



M^n" - P 



Act 



(91) 



Eq — Ueq _ 
+ to-UpQ — 



1 d 
7T d In r ' 



(tVt^ - 1) 



2 r d / y/i^ - 1 

3 7 d In r I r 



6P0 



7 
U, 



cq 



eq 

4 1 d , 



PO - C/eq/3 6E0 16 r d / 



eq 



where 



to 

h 

t2 



3 A5 ^ dlnr \ r I 
7 ^ dkP(kPfi^{k^/T) 



r 


2neaT^ 


_7 


vr 


r 


2neaT4 


_7 


vr 







dP(P 


fi^{k^/T) 


CTa 


+ CFs 


dk^{k^ 


)^^P{kVT) 



r2neaT^Jo aa + crs{l - a/b) 



(92) 
(93) 
(94) 

(95) 
(96) 
(97) 
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with Ucq = 7r^T^/15 as before. 

From the analytical model as well as the numerical results we know that 7T = constant and 
7 oc r are excellent approximations for r < r-y. Then 6F0 = 0, 6E0 ~ 0(to/7^) and 6P0 ~ 0(t2/7^) 
when 7 1; since t^^ and are of the same order as the absorption and total optical depth, 
respectively, we see that the corrections to blackbody radiation are negligible when 7 1. 

These results are also consistent with a more careful treatment of the equation of radiative 
transfer. This equation is most conveniently written as 



d 2(1- A) /rd7 A ( d 
dr r V7 dr ) \ dq dX 



7A 



-Mq,X,r)+efBB{qhT) + {l-e) T dX' f^{q, X' ,r) 

Jo 



(98) 



in an extremely relativistic flow (Grimsrud 1998). Here e = Oajioa + o"s), r] = 27^(1 — ^u), 
A = 1/(1 + r/) and q = k/2X, with fi the direction cosine and k the photon energy in the lab frame 
Again, the contribution from pairs is neglected on the right hand side of the transfer equation. 
The left hand side of this equation simplifies considerably when 7 oc r. If, in addition, 7T = 
constant, we see that /^(g, A,r) = fBB{q/jT) is the solution for the distribution function. Thus, 
as long as these two conditions are satisfied, the deviation from blackbody radiation vanishes in 
the high-7 limit up to terms of order 0(1/7^) which we have neglected. Note that equation (98) 
holds for both large and small optical depths. 



It is interesting to contrast the high-7 limit of equations (|92|)-(|94D with a static atmosphere, 
where 7 = 1. In the static case, neglecting pair- heating, 

5^0 = 0, SFo = -Ui^^, and 6P0 = (99) 
3 dmr 

in the diffusion approximation, r ^ 1. Often one relates the radiation pressure to the radiation 
energy via a variable Eddington factor /Edd (''")) Po = Eq[1/3 + /Edd(T")]- This is a closure relation 
needed for the first two moment equations obtained from the equation of radiative transfer. The 
above relations constrain the form of /Edd (''")• For r » 1, /Edd(''") must be at least quadratic in 
1/r; the linear term vanishes in this limit. 



Recall from Section 3.4 that a free-streaming radiation field, coupled to the matter via the 



dynamical equations T'^.f.f^ = —G", resulted in 7 oc r and consequently a radiation field close to 
blackbody out to optical depths much less than 1. This behavior is then also consistent with the 
full radiative transfer equation as discussed above. Furthermore, note that since deviations from 
blackbody in the radiation fields start building up only when r <C 1, we cannot use a variable 
Eddington factor /Edd(''") in our extremely relativistic atmosphere. 
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5.2. Observed spectrum 

To compute the observed spectrum, we use the fact that the (comoving frame) photon 
distribution function is blackbody to a very good approximation at any r for 7 oc r. This allows 
us to find the radially directed flux per energy interval seen by a stationary observer at arbitrary 
r, using 

r (100) 

dp 27r^ J-i exp[7p(l — vfi)/T] — 1 

where 7 oc r (and p is the observer-frame photon energy) and T is the blackbody temperature 
in the frame moving radially with Lorentz factor 7. (The same results can be found from free 
streaming at r < 1.) For 7^1, the integral can be done analytically to give 

At small values oip/^T, this expression tends to 

dF _ p'^T 
dp 27r^7 

at large values of p/^T it tends to 

dF p'^T 



ln(27r/p); (102) 



dp 2tt'^j 



exp(-p/27r). (103) 



In order to compare ( |10lD with a blackbody that fits the observed flux best, we minimize the 
integral 

S{Bo, xo) = rdxl -x^ ln[l - exp(-x)] ^ \ (104) 

Jo [ exp(x/xo) - 1 J 

with respect to the two parameters Bq and xq. The two functions are plotted in Fig. ^. Compared 
to a blackbody, the observed flux is seen to have a broader peak, and a shallower slope at low 
photon energies. A similar result was found by Goodman (1986) from his fireball simulations. 



5.3. Perturbations about spherical symmetry 

So far we have assumed the flow to be static and spherically symmetric. In reality, however, 
these assumptions break down if the flow is affected by, e.g., a magnetic field, an anisotropic 
and/or variable temperature distribution at the inner surface r = or inhomogeneities in the 
external medium. Here we will merely look at the effects on a static flow due to perturbations 
about its spherically symmetric solution. 

For simplicity, we assume the energy density and pressure to be functions of temperature 
only, i.e. we use the equilibrium expression for the number density of electrons and positrons. As 
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we have seen, this is a good approximation inside the photosphere as far as the dynamics of the 
flow is concerned. Also, the flow is assumed to be azimuthally symmetric. 

Defining the quantity a = {p + P) /T, the energy-momentum tensor can be written as 

T''^ = TaU'^U'^ + Pg'^^. (105) 
Contracting the equations of motion (T^^ = 0) with Ua results in the conservation law 

{aU^.^^ = 0. (106) 
This can be used to simplify the equations of motion; 

(g'^P + [/"C/^) Tp + TU^U^p = 0. (107) 



We now introduce the four- vector = TU", which allows us to write equation (107) as 

{Va-p - Vp.a) = V'^UJo.p = 0. (108) 

We have here defined the vorticity tensor oj^p = {TUa).fj — {TUp).^ = {TUq) p — (TUp) ^ . It is 
straightforward to show that the circulation 

{TUa).0 - (TUfs)] dx°dx^ = [ Lo^pdx'^dx^ (109) 



is unchanged along the path of a fluid element. 

We are interested in axisymmetric {d/d(j) = 0) perturbations around spherically symmetric 
flow. Let the temperature of the unperturbed flow be T(r), and the nonzero components of the 
unperturbed four-velocity be U^{r) = 7(r) and TT {r) = 7(r)u(r); consequently {r) = T{r)U^{r) 
and V^(r) = T{r)lT~ {r), and oJaji = for the unperturbed flow. Denote perturbed quantities by 
(5r(r, 9), 5U^{r, 9) and 5y (r, 9). From the fluid equations V^ajf^u = we find, to first order in the 
perturbations, 

V^UJf,o + V"uj^r=0. (110) 

The // = and fi = r components of this equation both imply ujQr = dSVo/dr = 0, so that 
= 7(^)5 independent of r. Below, it will prove convenient to expand 

f{0) = —Y,l{i + i)nPm. (Ill) 

where Pt{9^ is the Legendre function. We omit £ = in the sum because it may be absorbed 



into the background solution, for which Vq = constant. The n = 9 component of equation ( 110 ) 
becomes 

dSVe dSVr V'df{9) {vy^ dP,{9) 
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after substituting for 6Vo. The fj, = (j) component of equation ( 110 ) implies d5V^/dr = 0, so that 
= g{e), or sin^ OU'^' = g{e)/T{r)r^. 



Next, we perturb the entropy equation, which may be written in the form 

1 d{r^smeQ{T)V>') 



sin / 



where Q{T) = a{T)/T. To first order we find 



v 



d_ 

dr 



+ 



0, 



1 d{sme6Vg) 



sin ( 



de 







(113) 



(114) 



where u{T) = d\nQ{T)/d\nT , and 5V'^ = g'^^dVe = r~^5Ve. 

The normalization condition — = V^^V^ imphes, to first order in the perturbations, 

T5T = -V^f{e)-T5V''- 

using this relationship to eliminate 6V^ = 5Vr in equations (p. 121) and (114) implies 



(where we have also multiplied by sm.9) and 

+ sin^y 



d{siYie5Ve) , ^2 -.^ nTyrd[vR{'r){^T /T)] _ 2r sin0/(0)rV dlnT 



do 



dr 



dlnr 



(115) 



(116) 



(117) 



where i'R{r) = ^{T) — (T/V^)^. Difi'erentiate equation ( [116| ) with respect to 6 and equation (^ 
with respect to r to find 



V_d_ 



^ Q 



dr 



,5T 



1 d 

sin 6 do 



d f5T 



sm I 



de\T 



(118) 



rTTO 



where (recall that dV /dr = 0) 



rT dlnT 
{Vf dlnr 



(119) 



When the background fiow is extremely relativistic and dominated by extremely relativistic 
particles, we may take = V ,Ttxl/r, and = 2 up to corrections ~ 7^^. With these 
substitutions, H{r) = 1 and equation (|11^) becomes 



2 d 



" 2 d 


(5T\ 


1 d - 


'Wr 


Kt). 


sin 6 do 



sm 



<1 

do 



5T_ 



Y,l{i + l)nP,{0), (120) 
£^0 
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which has the general solutions 

6T 



J2[Ti + CiCosih/-fir)) + sesmike/jir))]Peie), (121) 



where — + l)/2, and q and si are constants to be determined from boundary conditions. 
Notice that since 7(r) oc r, the temperature fluctuations become independent of r asymptotically: 

,1™ 5- = E(^^ + (122) 



From equations ( |l2lD and ( 116 ) we find 



•^^e = 2^ - + 7— sm(/i;£/7(r)) - — cos(/i:£/7(r)) + — r; (123) 



the additional solution proportional to l/sin0 is a (singular) potential flow with no associated 
temperature perturbation. Since SVe = r'^T6U^ , and Tr = constant, the non-radial velocity 
component 6v^ = r6U^ /j oc at large radii (ignoring the singular term in eg. jl23[). From 



equation (115) we find that 




= E ^^^(^ + 1) + ^^hcosih/^{r)) + sesinikeMr))] P,{e)- (124) 



here, we ignored a correction term oc /{y^)"^ compared with i{i + 1) in the term oc r^. To 
leading order in 1/7^, SV^ = . We note that in this limit ujqq = ujQr = df{9)/d0. In the 
solutions given above, only the terms proportional to ti represent flows with nonzero vorticity; the 
remaining terms are the general solution for perturbed potential flow. 

Although memory of the temperature fluctuations at the inner boundary is maintained far 
out in the flow, it is the value of = 7T that determines whether there are any observable 
consequences. As a result, if f{0) = 0, so that the flow is derivable from a potential, the radiation 
spectrum detected by a distant observer is the same as for the background spherical flow. If 
f{9) 7^ 0, so that the flow has nonzero vorticity, distant observers detect a superposition of 
quasi-thermal spectra of the kind derived in Section While these results have been derived for 
small axisymmetric perturbations about spherically symmetric flow, it is conceivable that when the 
perturbations become significant, the observed spectrum can appear substantially non-thermal. 



6. Baryon Contamination and 700 

How is the M ^ limit achieved for nonzero (baryon) mass ejection rate Ml This question 
raises three subsidiary ones: (1) For a given E, what is the critical value of E/ M above which the 
flow is essentially the same as for M = 0? (2) Since we know that the asymptotic Lorentz factor 
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7oo is finite for M = (equation |7T|), and that 700 ~ E/M for relatively large M (small E/M), 
does 7oo grow monotonically with increasing EjMl (3) What is the largest possible value of 700? 

Below, we demonstrate that 700 ^ E/M until it reaches a critical value 



4:) ^?,bQ{Z/Af'\n/lQ^cxaflSf\Ti/me), (125) 
c,M 



M 



at which the radius vm ~ {E/M)ri where the flow becomes matter dominated first moves out 
to the photosphere. As M decreases further, we find that 700 remains virtually constant at 
~ {E / M)c^M until a second critical value 



M 



~ 7 X 10*(Z/.4)(ri/10'^cm)(ri/me), (126) 

C,P 



above which the densities of electrons and positrons are nearly equal once pairs go out of 
equilibrium. Between (£'/M)c,p and 



M 



~ 5 X 10^(ri/10^cm)(Ti/me) (127) 

c,0 



the inertia in pairs beyond the photosphere becomes progressively more important compared 
with the inertia in baryons, and 700 grows from ~ [E / M)c^m to roughly its maximum value, the 
M ^ limit, which applies at E/M > (S/M)c,o- Over the entire domain E/M > (E/M)c,m, the 
value of 7oo only grows by a factor ~ {Amp/2Zme)^/^ ~ 10. 



6.1. Equilibrium flow 

First, we review the results for an equilibrium flow with baryons (Paczyhski 1986). From 
conservation of baryon number and energy, it follows that 

N = Anr'^n'jv = M/Arup = constant; (128) 
E = 47rr2(p + P)72t; = constant. (129) 

As before, we get 7 oc r and T oc 1/r in a radiation dominated flow. This behavior is maintained 
until /3 + p ~ nArUp, which occurs at a radius r^/ ~ {E/M){ri/ji) where 7 ~ E/M. In an 
equilibrium flow at sufficiently low T the optical depth is essentially r ~ ZnaTf/'j, resulting in 

^ ^ ZN'^jAmpf-fjaT / {ru/rf for rM > r ^^^^^ 
inriE^ \ (rM/r) for vm < r. 

A critical value of E/M is defined as the value which gives r = 1 at r = vm, i-e. rph = vm- With 
p + P (11/4) (4/3) (7r2/15)i;'^ at this gives 
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The optical depth can then be rewritten as 



E/M 



c,M 



E/M 



{ru/r) 



for rjv/ > r 
for rM < r. 



(132) 



If E/M < {E/ M)c,M, fM is inside the photosphere and we expect there to be httle acceleration 
outside rM, so the asymptotic Lorentz factor is E/M. If, on the other hand, E/M > {E/M)c^m, 

r • . • • -|4/3 

then vm is outside the photospheric radius ~ rM (E / M)c^m/{E / M) , and the Lorentz 

factor at rph is 7ph ~ 7M(?'ph/'^A/) ~ / ^)t^M(-^ / ■^)~^^^ ■ "^^^ acceleration of the flow by the 
radiation flux beyond rph is found in a similar way as for the non-equilibrium pair wind, using 



Am, 



d'-f 
' dr 



ZaTFo{r) 



3 V 



r'j 



1 



(133) 



For rM > '''ph this results in a final Lorentz factor 700 ~ {E/AI)cM < E/M. 



The use of a free streaming form for the radiation field for r < 1 implies that the radiative 
component of E is constant in this regime. This is a good approximation when E/M > {E/M)cfl 
and the energy content of the baryons is negligible. When E/M < {E/AI)c^m, the radiation energy 
is much smaller than the baryon energy outside the photosphere, and so the details of the radiation 
field are not dynamically important. In the intermediate regime, radiative acceleration increases 
the baryonic energy, causing a corresponding decrease in the radiative energy unaccounted for 
when using the free streaming approximation. The error thus made can be estimated as 

1/3^ 

Ai?baryon/i? = (7oo " lpl.)M/E ~ ^^llJ^Ujb^ h 



E/M 



{E/M) 



c,M 



E/M 



(134) 



The maximum A£'baryon/-E' is ~ 0.1 and occurs at E/M = (4/3)^(£^/M)c,m- Numerically we 
find the maximum gain in baryonic energy to be about 20 per cent of the total flow energy, with 
a functional dependence on E/M following roughly the estimate in equation ( [1341) . Since the 
error made in assuming a free streaming radiation field is relatively small as far as the dynamics 
is concerned and only affects a small interval in E/M, we adopted this simplification in our 
numerical calculations. The resulting values of 700 are very accurate for E/M well above and well 
below (^/M)c,M, and should be correct to 10 - 20 per cent for E/M ^ {E/M)c,m- 



6.2. Baryons and non-equilibrium effects 



With baryons in the flow, charge neutrality requires that the electron number density n_, the 
positron number density rz-|_, and the baryon number density n obey n_ = n-|- + Zn. The positron 
number density n+ is then found from 

(^^^^+7^^) = -(o-annt') (n+ + Zn) n+ 



n, 



eq 



(135) 
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where, as before, 



"cq ^ (^+"'-)eq- At large M, the excess of electrons over positrons is always 
large at T < mg. The positron density first becomes important at rgq when 



E Znie jiTi C 



i 3 X 10'^—-—^ f ] (136) 

M Arup rUe (InC)'^ ArUe \W^cmJ 

where C = ■K{e^ /m1)m^ri/^i ~ 4.4 x 10^^(rj/10®cm)7~^; positrons first become important at the 
photosphere and beyond when 



E 
M 



E 
M 



C,P 



Zme 7»^i C 
Anip nif, (InC)" 



For E/M > (£'/M)c,p, the inertial mass density outside the photosphere is 

p ~ AmpTi + 2men+ = n+meff 
where the effective mass is approximately 



l + {M/E){Ti/m,)C^,/A{\nCf 



which varies between ~ Amp/Z for E/M ~ (£'/M)c,p and 2me for 



E 
M 



> 



E_ 

M 



c,0 



4me {In Cy- 



(137) 

(138) 
(139) 

(140) 



Notice that {E/M)c,o ~ (Amp/2Zme)(^/M)c,p. 

The acceleration caused by the radiation force for r < 1 can be estimated using 



d'j Att f rph7ph \ 



^ph7 



resulting in 



7o 



{Amp/2Zmef'\E/M),^ 



M 



1/4- 



(141) 



(142) 



1 + {M/E){Ti/me)C-ii/A{\nCY 

Equation ( |142D shows that 700 rises slowly (~ [E/MY/^) from ~ (^/M)c,m at E/M ~ (£^/M)c,p, 
to ~ {Amp/2ZmeY/'^{E/M)c,M at E/M > {E/M)c,o- Since 700 ~ (-E/M)c,m for 
{E/M)c,p > E/M > {E/M)c,M, the total increase in 700 for all E/M > {E/M)c,m is a 
factor ~ {Amp/2Zme)^^^ ~ 6. This alleviates any 'fine-tuning' problem necessary to produce large 
asymptotic Lorentz factors (although proportionality to Ti/rUe remains): Values within an order 
of magnitude of one another are found as long as E/M is sufficiently large. 

This qualitative behavior is seen in the numerical solution: 700 is shown as a function of M /E 
in Fig. clearly showing the four distinct regions discussed above. The transitions between the 
regions are seen to correspond to the critical points {E / M)c^m-, iE/M)c,p, and {E/AI)cfl, shown 
as dotted vertical lines in Fig. 0. 
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7. Discussion 

In order to make our treatment of the relativistic e^7 wind somewhat more complete, we 
address briefly (and largely qualitatively) two additional questions: (1) What happens if the initial 
temperature is much higher than the electron mass, e.g. comparable to the mass of muons or even 
nucleons? (2) Do weak magnetic fields alter the flow significantly? 



7.1. The effects of Ti > m-e 

When T ~ m^, the muons are in equilibrium with the radiation, + ^u" ^ 7 + 7. As T 
drops below a few percent of m^, the muon pair annihilations freeze out in the same way as the 
electron-positron pair annihilations do at T ~ O.OSmg. Going through the same calculations as 
for the electrons, we find that req,e/?'eq,^ ~ "i^/m-g, and that k, for r > req,e assuming that 
the muons have not decayed before reaching ^eq,^. (For req,^ <^ r <^ ''eq.e, n-e T <C and 
therefore <^ pe in this regime.) However, the mean lifetime of a muon is only ~ 2.2 x 10^^ 
seconds. They will therefore be significantly abundant at rgq^g only if > J^P'" dr/j{r), i.e., if 
7j > 80(rj/10^cm). For 7^ smaller than this, the presence of muons in the fiow will not have any 
significant infiuence on its dynamics. On the other hand, if 7^ is large enough for the muons to 
survive until 7 reaches 700, then they will reduce the asymptotic Lorentz factor slightly because 
of their added contribution to the inertia of the fiow. (Demanding > /jT^'*^ dr/7(r) requires 
7i > 1.9 X 10^[ri/10^cm].) For T <^ rUe, equation (p6|) reduces to 



(pe + Pf^)^rr^ « « -2ne{aa + as)Fo (143) 

amr 7 7 

when muons are abundant. (Since (Tt oc m~^, we may neglect the muon contribution to the 
scattering cross section.) Thus the radiation force is approximately equal to that of a fiow with 
only electrons and positrons, whereas the inertia is doubled since + ~ 2pe ~ 2mene in 



the region outside the photosphere. In the notation of Section 3^, then, A A/2. And since 
7oo ~ 7phA^''^, the muons will reduce the asymptotic Lorentz factor by a factor 2^/^ ~ 1.2, 
provided that most of them have not decayed by the end of the acceleration epoch. Ultimately, 
the relativistic muons decay, resulting in an additional non-thermal population of electrons. 

In the arguments above we implicitly assumed that muons and electrons are coupled. Since 
their coupling via photons is very weak, the relevant remaining mechanism is through Coulomb 
scattering. The time-scale for Coulomb interaction is roughly (See e.g. Spitzer 1978) 




m, 



2^2 



(3/2e^)Vr3/7rn, 



(144) 



Comparing this to the expansion time-scale, texp ~ r/^, we get tc/^exp ~ 7.2 x 10 ^(me/r)^/^, 
which is less than unity for T/nif, > 1.7 x 10"'^. Recall that Tph/me ~ 0.037 and T^/rrif, ~ 8 x 10^^; 
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this indicates that the electrons and muons are marginally coupled for most of the acceleration 
regime < r < r^. The reduction of 700 by a factor 1.2 therefore represents an estimate for the 
maximum change in the asymptotic Lorentz factor caused by the presence of muons in a flow with 
very high 7^ (or small rj). 

For high temperatures the inner portion of the wind can be optically thick to neutrinos as well. 
The cross section for electron-neutrino scattering is ~ cro(T/me)^, where ctq = 1.76 x 10~^^cm^. 
The corresponding optical depth is 

r^r) ~ / drn^ir)aAr)Mr) ~ 1.4 x lO^-^ (-) ; (145) 

Consequently, Tj^ > 1 for T/nie > 377j^^^(rj/10^cm)~-'^/^. Neutrinos are kept in thermal 
equilibrium via reactions like fe + i^e ^ e~ + e"^, whose rate exceeds the expansion rate for 
T/rrie > 3l7,^''^(ri/10'5cm)-V5. 

If T increases further we have to include not only muons in the flow, but also mesons. However, 
since mesons are very short-lived, their dynamical effect on the flow will be negligible. But for 
initial temperatures high enough for nucleon-anti-nucleon pairs to exist, there is the possibility 
that the surviving fraction (after annihilations freeze out) will result in a significant baryon loading 
in the wind. The cross section for nucleon-anti-nucleon annihilations is (av) ~ l/rn^, resulting in 
a freeze-out temperature of Teq nuci/'n'p ~ 0.026. This results in negligible baryon loading in the 
wind: Pnud/ Pe± ~ 4 x 10"^ for r > req,e- 



7.2. Effects of weak magnetic fields 

In this subsection we will estimate how weak magnetic fields can affect the flow. The 
'background' wind is assumed to expand radially, and we consider the effects of radial and 
tangential fields. 



In Section 5.3 we presented a general treatment of perturbations of a flow with equation of 
state that depends only on temperature. When magnetic fields are present, and the electric field 
vanishes in the rest frame of the flow, the equation of entropy conservation, {crU^)-f^ = remains 
true, but the remaining fluid equations are modified to 

where F^a is the Maxwell tensor, and J'^ is the current density four vector. 



We can use equation ( |146| ) to examine the effects on the flow of a magnetic field; to do so 
requires solving Maxwell's equations simultaneously. A particularly simple example is a radially 
directed magnetic field, as might arise, for example, if the flow originates on a magnetized star, 
and pulls magnetic field lines outward along with it. In that case, it is easy to show that, assuming 
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axisymmetry, the magnetic field strength is proportional to b{9)/r'^, where b{9) is arbitrary. To 
lowest order, the electric field associated with a radial magnetic field vanishes in the frame of a 
stationary observer as well as in the frame comoving with the flow. 

Under such conditions, the magnetic field engenders no change in 7T to first order. This 
follows from the n = r component of equation (|146|), which implies V^wor = 0, or dSVo/dr = 0. (A 



similar conclusion follows from the = component of eq. jl46[ .) Thus, the only perturbations to 
6V0 are those imposed at the inner radius of the flow, as discussed in Section ^.3| none are created 
by the field. Radial magnetic fields drive non-radial flows whose amplitudes decline with radius 
even outside the photosphere (Grimsrud 1998). 

Changes in 7T require a non-radial magnetic field. The condition of perfect conductivity, 
E -|- V X B = 0, and Faraday's law imply E = —'Vip = — v x B in steady state; for axisymmetry 
we find, to lowest order, 

^ = rr ^ = y^'^') 

rv(r) at) r at) 

Associated with ordered tangential fields is a Poynting flux in the stationary frame, 



(148) 



there is no Poynting flux for radially directed B to lowest order (since E = 0), and there is, of 
course, no Poynting flux in the comoving frame. Note that Ai^r'^S decreases with increasing radius 
as v{r) — > 1 from below. For these fields the /i = component of equation (|146| ) implies 



dWo _ 1 f dipjeW d n\ 1 / #(g) y d7 



dr 47rr2Q(r)y' V dO J dr \v J 47rr2Q(r)y V^H dO ) dr' ^^^^^ 



the n = r component of equation ( |146| ) yields an identical result. In the extreme relativistic limit, 
where 7 oc r and Q{T) (x T (x r^'^, equation (|149|) has the solution 



1 / #(g) y _ |B| 

87rr2Q(T)7217° V dO ) ~ 87rQ{T)^W°' 



so 6Vq oc r^2_ Although tangential fields alter 7T, the perturbation peaks near the lift-off radius, 
and decreases far out in the flow. We therefore expect little or no effect on the observed radiation 
spectrum as a consequence of such fields. There is additional radial acceleration of the pairs 
outside the photosphere as a consequence of the tangential field; we estimate the change in bulk 
Lorentz factor of the pairs to be A7 « {dil>{6) / dO)'^ / AnieN'^'^^ if pair annihilation can be neglected 
outside the photosphere, so the pair loss rate N becomes independent of radius. 
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8. Conclusions 

Radiation energy can escape from a fireball in two different ways: If there is significant baryon 
contamination present, much of the energy will be converted into bulk kinetic energy (Shemi & 
Piran 1990). However, when the expanding atmosphere has swept up a significant amount of 
surrounding matter, kinetic energy can be converted into escaping radiation at the resulting shock 
front (Meszaros & Rees 1993). In a similar way, internal shocks due to a non-uniform velocity can 
convert kinetic energy into radiation (Rees & Meszaros 1994). The other mechanism for radiation 
escape is more direct: If the particle content is small, the fireball can become optically thin before 
being matter dominated. 

In this paper we have considered an extreme case of the latter possibility, in which there are no 
baryons present. The opacity is then due to electron-positron pairs, resulting in a very large optical 
depth for temperatures greater than the electron mass. Further out in the flow the temperature 
decreases, pair creation is suppressed and annihilations freeze in; this results in a small but 
non- negligible amount of surviving pairs. The radiation force acting on the particles accelerates 
the pairs considerably, even after the flow has become optically thin. We found that the Lorentz 
factor of the flow approaches the constant value 700 ~ 

only when the optical depth falls below ~ 1.7 x IQ^^^fJ^ . This increases the asymptotic 
energy content of the pairs by a large factor, their fraction of the total energy approaching 
~ 8.5 X lQ^^^l^^[{(Ta/ (^t) + {'^s/ (^t)]^^^ ■ The flow is always radiation dominated for reasonable 
values of the input parameters. 

Even if the initial temperature were much higher than the electron mass, the resulting flow 
would not deviate significantly from an e^7 wind. If there are muons present, they will decay 
before the electron-positron annihilations freeze out, unless 7^ is very large. And even for 7j high 
enough for the muons to survive until far outside the photosphere, their added inertia will only 
reduce the asymptotic Lorentz factor by at most 20 per cent. For even higher Tj, one may have 
nucleon-anti-nucleon pairs present in the flow. However, nucleon-anti-nucleon annihilations freeze 
out at a relatively low temperature, thus causing the baryon contamination in the flow to be 
negligible. 

The photon distribution function in the comoving frame is very close to that of blackbody 
radiation. This is because 7 oc r and 7T ~ constant are excellent approximations in the flow 
until r = Vy where the optical depth is t-j ~ 1.7 x 10~^7j^^^. Practically all the observed radiation 
therefore originates from a region where these two approximations hold. As was discussed in 
Section |5|, the conditions 7 tx r and 7T = constant imply that the equation of radiative transfer is 
solved by a blackbody distribution function in the comoving frame of the flow. The spectrum seen 
by an observer in the lab frame will deviate somewhat from a blackbody in that it has a broader 
peak and a shallower slope at low photon energies. Such spectra are not typical of observed 
7— ray bursts, which are characterized by flat fluxes for logarithmic energy intervals (e.g. Schaefer 
et al. 1992, 1994; Kouveliotou 1994). A superposition of quasi-thermal spectra from numerous 
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source regions radiating independently (but with different physical parameters) might produce flat 



spectra. The non-radial perturbation calculations discussed in Section |5^ lend partial support to 
this idea. 

Ordered magnetic fields (whose energy content is small compared to that of the flow) will 
make the temperature and the velocity of the flow anisotropic and may enhance the bulk Lorentz 
factor of pairs beyond the photosphere, but do not affect the spectrum seen by a distant 
observer significantly. 

The results obtained for zero baryon number also apply when the baryon loading is sufficiently 
small. For very large baryon loading, the flow becomes matter dominated at optical depths larger 
than one, and in this case the asymptotic Lorentz factor 700 ~ E/M. As E/M increases above 
{E/M)c^M ~ 350(.Z^/yl)^/^(rj/10^cm)^/^7j^''^rj/me, the asymptotic Lorentz factor at first levels off 
at 7oo ~ (£^/M)c,M. For stih larger E/M> {E/M)c,p ~ 7 x W^{Z/A){ri/W^cm){Ti/me), the 
asymptotic Lorentz factor rises ~ {E/Uy/^, until E/M ~ {E/M)c,o ~ iAmp/2Zme)iE/M)c,p. 
The M ^ limit applies for E/M > {E/M)c,o ~ 5 x W'^ {ri/ cm) (Ti/ me); in this regime, 
7oo ~ {Amp/2Zme)^^'^{E/ M)c^M- The fraction of the total wind luminosity that emerges in the 
form of bulk kinetic energy falls below one at E/M ~ (£'/M)c,M) and decreases monotonically 
until asymptoting to a finite value ~ 10~^7j^'^^ as M ^ (see equation [|7^). Thus, for all 
E/M > (E / M)c^M, the asymptotic Lorentz factor varies by a factor of only ~ {Amp/2Zme)^/^ ~ 6. 
The maximum possible 700 is the value found for M = and is finite. Although we have only 
considered steady winds here, it seems likely that similar results would hold for fireballs originated 
impulsively. 

This research was supported in part by NSF grant AST 93-15375 and NASA grant NAG 
5-3097. OMG thanks Nansenfondet for financial support. 
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Fig. 1. — Summary of analytical model with 7j and rj/lO^cm both of order unity. 
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Fig. 2. — The ratios rrajfi and rcc^/vm are shown as functions of the initial temperature Tj for fixed 
initial Lorentz factor 7^ in the upper panels, and as a function of 7^ for fixed Tj in the lower panels. 
Vi is the initial radius, = r{T = me) and rgq = r(ne = 2ne,eq)- The absorption cross section 
is (7a = 10~^crT- The solid lines correspond to the numerical solution, and the dotted lines to the 
analytical model. 
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Fig. 3. — The ratios rph/req and r-y/rph are shown as functions of initial temperature Tj and initial 
Lorentz factor 7^. r^q = r(ne = 2ne,eq)) ''ph = = 1) and r-y = r{( = 1/2). The absorption 
cross-section is fixed at cja = 10~^cjt- The solid lines correspond to the numerical solution, and the 
dotted lines to the analytical model. 
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10° IQi 102 IQ3 



T,/m^ 

Fig. 4. — The asymptotic Lorcntz factor as a function of the initial temperature for 7^ = 2 and 
7i = 10. The absorption cross-section is held fixed at aa = IO^^ot- 
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Fig. 5. — The local temperature T{r) in units of nie, and its derivative d In T/d In r(r), are plotted 
for different absorption cross-sections. In these plots, Ti = lOmg and 7^ = 2. 
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Fig. 6. — The asymptotic pair number density ng times = {r/rif' is shown in the left panels in 
units of rrZg as a function of initial temperature Tj and Lorcntz factor 7,. The right panels show the 
asymptotic energy content in the pairs relative to that in the radiation. In all plots, Ua = 10~^(7r- 
Solid lines represent the numerical solution and dotted lines the analytical model. 
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Fig. 7. — The Lorentz factor j{r) and its derivative dlnj/d\nr(T) are shown for a set of different 
absorption cross-sections. Tj = lOmg and 7i = 2 are fixed. 
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Fig. 8. — The spectrum seen by an observer in the lab frame (soHd Hne) compared to blackbody 
(dotted hne). The units for the flux are arbitrary. 




10-' 10-2 10-3 10-* 10-=> 10-8 10-' 10-8 10-3 10-1" 

ii/E 



Fig. 9. — Asymptotic Lorentz factor as a function of M/E, the ratio of the mass injection 
rate to the energy injection rate. The dotted vertical lines correspond to the transition points 
(^/M)c,M, (^/M)c,P, and (^/M)c,o. In this example, T^/me = l,7i = 2,a„/aT = IQ-^a^/aT = 1 
and A = Z = 1. 



